This project seeks to estimate sport fish harvest and releases of rockfish in Alaska waters by improving on the Howard et al. (2020) methods and expand the time series back to 1977 when the statewide harvest survey (SWHS) was first implemented. This is essentially a Bayesian version of the Howard methods that allows for more appropriate and defensible sharing of information between areas, handles missing data in a more appropriate manor, accurately propagates uncertainty throughout the estimation procedure and replaces the Howard decision tree approach to low sample sizes with a hierarchical model. The methods and results for generating harvest estimates are generally consistent between the Bayesian model and the Howard methods. Harvest estimates are consistent with Howard estimates during contemporary times, but may differ based on more appropriate weighting of SWHS and logbook data, including estimating and correcting bias in the SWHS data.

Table 1. Summary of key improvements in reconstructiing sport fish removals of rockfish using the Bayesian model as compared to the Howard et al. (2020) methods.
Issue Howard Bayes
Time series 1999 - present 1977 - present
Bias in SWHS Not explicitly dealt with. Relies on logbook data and ratios of guided/unguided from SWHS data to estimate unguided releases and harvests. Explicitly estimates bias in SWHS harvest and release estimates based on logbook data.
Species composition of releases Assumes that species composition of releases is equal to that of the harvest, which is not evident in the logbook data. Recognizes different release probabilities by species / species assemblage and estimates it from logbook data and bias corrected SWHS data
Sample size limitations Uses sample size threshholds such that when areas fall below those threshholds values are borrowed from nearby areas. Uses a hierarchichacal modelling approach that shares information between areas in the same region. Thus all data is used, even with small sample sizes. This is a more sound method that avoids assumptions and uses all of the data.
Error propogation Error is propogated when variance estimates are available, but there is uncertainty associated with borrowing values from nearby areas, or the assumption of species compositions being identical in harvest and releases, are not dealt with. By breaking the assumption that species composition is equal between harvests and releases, uncertainty in the release estimates is more reflective of the fishery. Furthermore, the hyerarchichal approach more accurately captures uncertainy within and between areas within a region.
**Figure 1.**- Data sources for estimating rockfish harvests and releases in ADF&G commercial fisheries management units. Note that initial rockfish harvest estimates were not differentiated into species assemblage or species until 1998 when logbooks began differentiating by pelagic and non-pelagic. Logbooks began to collect data on yelloweye beginning in 2006. Port sampling programs to gather data on species composition of harvests began in 1996 in Southcentral and Kodiak and in 2006 in Southeast.

Figure 1.- Data sources for estimating rockfish harvests and releases in ADF&G commercial fisheries management units. Note that initial rockfish harvest estimates were not differentiated into species assemblage or species until 1998 when logbooks began differentiating by pelagic and non-pelagic. Logbooks began to collect data on yelloweye beginning in 2006. Port sampling programs to gather data on species composition of harvests began in 1996 in Southcentral and Kodiak and in 2006 in Southeast.


**Figure 2.**- SWHS harvest (left) and release (right) estimates from guided trips (x-axis) versus repoted harvests from charter logbooks (y-axis).

Figure 2.- SWHS harvest (left) and release (right) estimates from guided trips (x-axis) versus repoted harvests from charter logbooks (y-axis).


Results

**Figure X.**- Rhat values and proportion of parameters that converged (Rhat < 1.1.)

Figure X.- Rhat values and proportion of parameters that converged (Rhat < 1.1.)

Estimate comparison

Since previous estimates of rockfish harvest have been produced these first 3 graphs will be used to show how the modeled estimates compare to the estimates produced earlier. For total rockfish the estimates are in general agreement although differences are noted. These estimates should be more reliable because they include both SWHS and guide logbook data, handle variance more appropriately, use hierarchical distributions when data is missing, directly consider observation error and are produced using reproducible research.

**Figure 3.**- Total rockfish harvests 1996-2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.

Figure 3.- Total rockfish harvests 1996-2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.



**Figure 3.**- Total rockfish releases 1996--2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.

Figure 3.- Total rockfish releases 1996–2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.


Notes from Adam: When looking at only black rockfish the most significant differences are for the Prince William Sound Inside area. I did not spend a great deal of time tracking this down although it looks like the previous version used bad values for \(P_{(black)ayu}\) for at least unguided anglers. For the moment I would ignore the results for BSIA and SOKO2SAP. I think it is possible to give approximate values for these areas but it will require a little more coding which I have yet to do.

**Figure 4.**- Black rockfish harvests 1996--2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.

Figure 4.- Black rockfish harvests 1996–2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.


And black rockfish releases…

**Figure 5.**- Black rockfish releases 1996--2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.

Figure 5.- Black rockfish releases 1996–2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.





**Figure 6.**- Yellow rockfish harvests 1996--2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.

Figure 6.- Yellow rockfish harvests 1996–2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.



**Figure 7.**- Yellow rockfish releases 1996--2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.

Figure 7.- Yellow rockfish releases 1996–2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.





**Figure 8.**- DSR rockfish (excluding yelloweye) harvests 1996--2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.

Figure 8.- DSR rockfish (excluding yelloweye) harvests 1996–2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.



**Figure 9.**- DSR rockfish releases (including yelloweye) 1996--2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.

Figure 9.- DSR rockfish releases (including yelloweye) 1996–2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.



**Figure 11.**- Slope rockfish harvests 1996-2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.

Figure 11.- Slope rockfish harvests 1996-2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.



**Figure 12.**- Slope rockfish releases 1996--2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.

Figure 12.- Slope rockfish releases 1996–2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.

**Figure 13.**- *Other* rockfish (non-pelagic and non-yelloweye) harvests 1996--2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.

Figure 13.- Other rockfish (non-pelagic and non-yelloweye) harvests 1996–2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.



**Figure 14.**- *Other* rockfish (non-pelagic and non-yelloweye) releases 1996--2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.

Figure 14.- Other rockfish (non-pelagic and non-yelloweye) releases 1996–2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.



Total Biomass Removal Estimates

**Figure 15.**- Black rockfish estimated total removals in lbs in 1996--2023. Total removals are calculated from the harvests, the estimated release mortality and mean fish weights. Error polygons represent 95% confidence intervals.

Figure 15.- Black rockfish estimated total removals in lbs in 1996–2023. Total removals are calculated from the harvests, the estimated release mortality and mean fish weights. Error polygons represent 95% confidence intervals.



**Figure 16.**- Yellow rockfish estimated total removals in lbs in 1996--2023. Total removals are calculated from the harvests, the estimated release mortality and mean fish weights. Error polygons represent 95% confidence intervals.

Figure 16.- Yellow rockfish estimated total removals in lbs in 1996–2023. Total removals are calculated from the harvests, the estimated release mortality and mean fish weights. Error polygons represent 95% confidence intervals.

**Figure 17.**- Pelagic rockfish estimated total removals in lbs in 1996--2023. Total removals are calculated from the harvests, the estimated release mortality and mean fish weights. Error polygons represent 95% confidence intervals.

Figure 17.- Pelagic rockfish estimated total removals in lbs in 1996–2023. Total removals are calculated from the harvests, the estimated release mortality and mean fish weights. Error polygons represent 95% confidence intervals.


**Figure 18.**- Non-yelloweye, demersal shelf rockfish estimated total removals in lbs in 1996-2023. Total removals are calculated from the harvests, the estimated release mortality and mean fish weights. Error polygons represent 95% confidence intervals.

Figure 18.- Non-yelloweye, demersal shelf rockfish estimated total removals in lbs in 1996-2023. Total removals are calculated from the harvests, the estimated release mortality and mean fish weights. Error polygons represent 95% confidence intervals.


**Figure 19.**- Slope rockfish estimated total removals in lbs in 1996-2023. Total removals are calculated from the harvests, the estimated release mortality and mean fish weights. Error polygons represent 95% confidence intervals.

Figure 19.- Slope rockfish estimated total removals in lbs in 1996-2023. Total removals are calculated from the harvests, the estimated release mortality and mean fish weights. Error polygons represent 95% confidence intervals.


Model fit

Logbook residuals

**Figure 20.**- Residuals from logbook harvests.

Figure 20.- Residuals from logbook harvests.


SWHS residuals

**Figure 21.**- Residuals from SWHS harvests.

Figure 21.- Residuals from SWHS harvests.



**Figure 22.**- Residual of SWHS releases.

Figure 22.- Residual of SWHS releases.

Parameter estimates

P(Charter)

These histograms show the posterior distribution of the mean percent of rockfish harvested by the charter fleet.

**Figure 23.**- Mean percent of harvest by charter anglers.

Figure 23.- Mean percent of harvest by charter anglers.


When considered annually we see the percent of rockfish harvested by the charter fleet follows our data fairly well although the model smooths out the changes and we just do not have much information about this ratio. Prior to 2011 the percent charter is confounded with SWHS bias and should be mostly discounted.

**Figure 24.**- Annual estimates of the percent of harvest by charter anglers for 16 commerical fishing manamgent areas, 1996-2023.

Figure 24.- Annual estimates of the percent of harvest by charter anglers for 16 commerical fishing manamgent areas, 1996-2023.

P(Harvest)

These plots show the fitted logistic line to the proportion of caught rockfish that are harvested. These estimates are used for hindcasting catch estimates based on the harvest data in early years when catch estimates are unavailable.


**Figure 25.**- Annual proportion of pelagic rockfish catch that was harvested. Note that pre-1990 estimates are used to estimate catch in these years when catch estimates are not available.

Figure 25.- Annual proportion of pelagic rockfish catch that was harvested. Note that pre-1990 estimates are used to estimate catch in these years when catch estimates are not available.


**Figure 26.**- Annual proportion of yelloweye rockfish catch that was harvested. Note that pre-1990 estimates are used to estimate catch in these years when catch estimates are not available.

Figure 26.- Annual proportion of yelloweye rockfish catch that was harvested. Note that pre-1990 estimates are used to estimate catch in these years when catch estimates are not available.


**Figure 27.**- Annual proportion of non-pelagic, non-yelloweye rockfish catch that was harvested. Note that pre-1990 estimates are used to estimate catch in these years when catch estimates are not available. Note, that this is not estimated for Southeast areas because non=pelagics are divided between DSR (including yelloweye) and Slope species.

Figure 27.- Annual proportion of non-pelagic, non-yelloweye rockfish catch that was harvested. Note that pre-1990 estimates are used to estimate catch in these years when catch estimates are not available. Note, that this is not estimated for Southeast areas because non=pelagics are divided between DSR (including yelloweye) and Slope species.


**Figure 28.**- Annual proportion of DSR rockfish catch that was harvested. Note that pre-1990 estimates are used to estimate catch in these years when catch estimates are not available. Note that the observed logbook data is for all non-pelagic, non-yelloweye fish.

Figure 28.- Annual proportion of DSR rockfish catch that was harvested. Note that pre-1990 estimates are used to estimate catch in these years when catch estimates are not available. Note that the observed logbook data is for all non-pelagic, non-yelloweye fish.


**Figure 29.**- Annual proportion of Slope rockfish catch that was harvested. Note that pre-1990 estimates are used to estimate catch in these years when catch estimates are not available. Note that the observed logbook data is for all non-pelagic, non-yelloweye fish.

Figure 29.- Annual proportion of Slope rockfish catch that was harvested. Note that pre-1990 estimates are used to estimate catch in these years when catch estimates are not available. Note that the observed logbook data is for all non-pelagic, non-yelloweye fish.

SWHS bias

Figure 23 shows the mean estimate for SWHS bias in harvests and releases. Cook Inlet, North Gulf Coast and North Southeast Inside all look pretty good while most other areas have substantial bias. Prince William Sound Inside has the largest bias. Bias in release estimates is substantial and whereas the SWHS appears to underestimate harvests, it appears to greatly overestimates releases by a factor of 2 or more in most areas as derived from logbook reported releases.

**Figure 30.**- Mean SWHS bias for harvests and catches. Note that a bias < 1 indicates that the SWHS *underestimates* the true value and bias > 1 indicates the survey *overestimates* the true value.

Figure 30.- Mean SWHS bias for harvests and catches. Note that a bias < 1 indicates that the SWHS underestimates the true value and bias > 1 indicates the survey overestimates the true value.


Our estimates of SWHS harvest bias track observations fairly well when he have guided harvest estimates. The estimates of release bias in the SWHS data track observed patterns to an extent, but appear to smooth these more volatile disagreements with the logbook data. Adam postulated in his initial start on this that some of this could be the result of the estimates of the proportion guided. This value was not modelled with a trend and thus applies a constant estimate when hindcasting. Data on these relationships could greatly improve this model.

**Figure 31.**- Annual estimates of SWHS bias in harvests and releases for 16 commerical fishing manamgent areas, 1996-2023. Note that a bias < 1 indicates that the SWHS *underestimates* the true value and bias > 1 indicates the survey *overestimates* the true value.

Figure 31.- Annual estimates of SWHS bias in harvests and releases for 16 commerical fishing manamgent areas, 1996-2023. Note that a bias < 1 indicates that the SWHS underestimates the true value and bias > 1 indicates the survey overestimates the true value.

P(pelagic)

We model the percentage of pelagic rockfish in the harvest because we have the information for charter anglers (via logbooks) starting in 1998. Other than looking at the model estimates you can use Figure 25 to compare the two data streams for pelagic rockfish harvest. In general they are in agreement with major exceptions in Price William Sound inside, Prince William Sound outside (early in the time series) and South Southeast inside.

**Figure 32.**- Annual estimates of the percent of the sport harvest that was pelagic rockfish for 16 commerical fishing manamgent areas, 1996-2023.

Figure 32.- Annual estimates of the percent of the sport harvest that was pelagic rockfish for 16 commerical fishing manamgent areas, 1996-2023.

P(black|pelagic)

Note that in Southeast Alaska we only have composition data starting in 2006. Tania dug up old SE data, but it did not provide any useful data for species apportionment. For the most part, P(black|pelagic) is relatively constant across areas, with the exception of Cook Inlet and NSEI in Southeast AK. It may be worth discussing whether the shifts in those areas is a result of improved or changing species identification rather than actual shift in the species composition of the catch.

**Figure 33.**- Annual estimates of the percent of the sport harvest of pelagic rockfish that were black rockfish for 16 commerical fishing manamgent areas, 1996-2023. Kodiak panels include data from a hydroacoustic survey and the proportion of pelagic rockfish that are black in those areas (red) and the adjusted proportions based on obseved harvests for charter (blue) and private (cyan) users.

Figure 33.- Annual estimates of the percent of the sport harvest of pelagic rockfish that were black rockfish for 16 commerical fishing manamgent areas, 1996-2023. Kodiak panels include data from a hydroacoustic survey and the proportion of pelagic rockfish that are black in those areas (red) and the adjusted proportions based on obseved harvests for charter (blue) and private (cyan) users.

P(yelloweye|non-pelagic / yelloweye|DSR)

**Figure 34.**- Annual estimates of the percent of the sport harvest of non-pelagic rockfish that were yelloweye rockfish for 16 commerical fishing manamgent areas, 1996-2023. Note that P(yelloweye) is the the proportion relative to non-pelagics for Central and Kodiak areas but is relative to DSR for Southeast areas.

Figure 34.- Annual estimates of the percent of the sport harvest of non-pelagic rockfish that were yelloweye rockfish for 16 commerical fishing manamgent areas, 1996-2023. Note that P(yelloweye) is the the proportion relative to non-pelagics for Central and Kodiak areas but is relative to DSR for Southeast areas.

P(DSR|non-pelagic)

**Figure 35.**- Annual estimates of the percent of the sport harvest of non-pelagic rockfish that were DSR rockfish for 6 Southeast commerical fishing manamgent areas, 1996-2023.

Figure 35.- Annual estimates of the percent of the sport harvest of non-pelagic rockfish that were DSR rockfish for 6 Southeast commerical fishing manamgent areas, 1996-2023.

P(slope|non-pelagic)

**Figure 36.**- Annual estimates of the percent of the sport harvest of non-pelagic rockfish that were slope rockfish for 6 southeast commerical fishing manamgent areas, 1996-2023. Note that P(yelloweye) is the the proportion relative to non-pelagics for Central and Kodiak areas but is relative to DSR for Southeast areas.

Figure 36.- Annual estimates of the percent of the sport harvest of non-pelagic rockfish that were slope rockfish for 6 southeast commerical fishing manamgent areas, 1996-2023. Note that P(yelloweye) is the the proportion relative to non-pelagics for Central and Kodiak areas but is relative to DSR for Southeast areas.



P(slope|non-pelagic & non-yellowye) For release estimates

## NULL



Weight Fits

**Figure 37.**- Mean weights of black rockfish used to estimate biomass from harvest and release mortality estimates in numbers of fish.

Figure 37.- Mean weights of black rockfish used to estimate biomass from harvest and release mortality estimates in numbers of fish.


**Figure 38.**- Mean weights of yelloweye rockfish used to estimate biomass from harvest and release mortality estimates in numbers of fish.

Figure 38.- Mean weights of yelloweye rockfish used to estimate biomass from harvest and release mortality estimates in numbers of fish.


**Figure 39.**- Mean weights of non-black, pelagic rockfish used to estimate biomass from harvest and release mortality estimates in numbers of fish.

Figure 39.- Mean weights of non-black, pelagic rockfish used to estimate biomass from harvest and release mortality estimates in numbers of fish.


**Figure 40.**- Mean weights of non-yelloweye, demersal shelf rockfish used to estimate biomass from harvest and release mortality estimates in numbers of fish.

Figure 40.- Mean weights of non-yelloweye, demersal shelf rockfish used to estimate biomass from harvest and release mortality estimates in numbers of fish.


**Figure 41.**- Mean weights of slope rockfish used to estimate biomass from harvest and release mortality estimates in numbers of fish.

Figure 41.- Mean weights of slope rockfish used to estimate biomass from harvest and release mortality estimates in numbers of fish.


Summary of unconverged parameters:

The proportion of unconverged species composition proportional estimates that are near the bounds of 0 or 1 is 0.5204082

The proportion of unconverged true parameters (defined as Rhat > 1.01) that are also in excess of an Rhat of 1.1 is 0.1 while 0.9 are less than 1.1.

Table Xa. Species proportion derived quantities that are not with 0.05 of the 0 and 1 bounds and are not converged relative to an Rhat of 1.01.

variable

area_n

year

user

mean

median

lower

upper

Rhat

sd

skew

area

p_pelagic

13

1

1

0.301

0.277

0.103

0.616

1.144

0.135

0.178

NSEI

p_pelagic

13

1

2

0.124

0.105

0.033

0.318

1.131

0.075

0.255

NSEI

p_pelagic

13

10

1

0.317

0.299

0.110

0.615

1.118

0.133

0.136

NSEI

p_pelagic

13

10

2

0.132

0.115

0.036

0.330

1.138

0.076

0.225

NSEI

p_pelagic

13

11

1

0.314

0.295

0.113

0.605

1.113

0.132

0.144

NSEI

p_pelagic

13

11

2

0.135

0.114

0.038

0.333

1.147

0.078

0.268

NSEI

p_pelagic

13

12

1

0.326

0.310

0.117

0.620

1.108

0.134

0.119

NSEI

p_pelagic

13

12

2

0.135

0.117

0.039

0.339

1.163

0.078

0.232

NSEI

p_pelagic

13

13

1

0.326

0.311

0.118

0.607

1.113

0.131

0.115

NSEI

p_pelagic

13

13

2

0.137

0.120

0.039

0.330

1.137

0.076

0.225

NSEI

p_pelagic

13

14

1

0.309

0.292

0.121

0.587

1.098

0.121

0.140

NSEI

p_pelagic

13

14

2

0.135

0.120

0.040

0.314

1.124

0.073

0.206

NSEI

p_pelagic

13

15

1

0.283

0.267

0.102

0.557

1.127

0.120

0.134

NSEI

p_pelagic

13

15

2

0.130

0.115

0.040

0.312

1.124

0.070

0.214

NSEI

p_pelagic

13

16

1

0.321

0.302

0.128

0.597

1.109

0.124

0.153

NSEI

p_pelagic

13

16

2

0.143

0.128

0.044

0.325

1.117

0.074

0.203

NSEI

p_pelagic

13

17

1

0.334

0.321

0.139

0.595

1.080

0.120

0.108

NSEI

p_pelagic

13

17

2

0.149

0.133

0.047

0.334

1.088

0.074

0.216

NSEI

p_pelagic

13

18

1

0.345

0.331

0.144

0.611

1.073

0.121

0.116

NSEI

p_pelagic

13

18

2

0.155

0.141

0.050

0.343

1.064

0.076

0.184

NSEI

p_pelagic

13

19

1

0.347

0.332

0.144

0.625

1.051

0.124

0.121

NSEI

p_pelagic

13

19

2

0.155

0.142

0.053

0.336

1.057

0.074

0.176

NSEI

p_pelagic

13

2

1

0.300

0.278

0.099

0.603

1.126

0.135

0.163

NSEI

p_pelagic

13

2

2

0.126

0.106

0.035

0.314

1.135

0.074

0.269

NSEI

p_pelagic

13

20

1

0.357

0.345

0.163

0.616

1.058

0.119

0.101

NSEI

p_pelagic

13

20

2

0.164

0.150

0.054

0.355

1.058

0.077

0.182

NSEI

p_pelagic

13

21

1

0.374

0.364

0.167

0.645

1.045

0.125

0.080

NSEI

p_pelagic

13

21

2

0.173

0.160

0.059

0.371

1.039

0.080

0.162

NSEI

p_pelagic

13

22

2

0.172

0.159

0.065

0.354

1.026

0.075

0.173

NSEI

p_pelagic

13

23

2

0.111

0.101

0.041

0.237

1.017

0.051

0.197

NSEI

p_pelagic

13

3

1

0.298

0.275

0.097

0.606

1.161

0.137

0.168

NSEI

p_pelagic

13

3

2

0.125

0.106

0.034

0.326

1.189

0.077

0.247

NSEI

p_pelagic

13

4

1

0.306

0.286

0.101

0.619

1.151

0.135

0.148

NSEI

p_pelagic

13

4

2

0.124

0.103

0.034

0.331

1.153

0.077

0.274

NSEI

p_pelagic

13

5

1

0.301

0.279

0.097

0.607

1.151

0.134

0.164

NSEI

p_pelagic

13

5

2

0.127

0.109

0.034

0.321

1.163

0.076

0.236

NSEI

p_pelagic

13

6

1

0.307

0.288

0.102

0.615

1.146

0.135

0.141

NSEI

p_pelagic

13

6

2

0.126

0.108

0.033

0.320

1.180

0.075

0.239

NSEI

p_pelagic

13

7

1

0.307

0.288

0.109

0.607

1.137

0.133

0.143

NSEI

p_pelagic

13

7

2

0.128

0.108

0.034

0.326

1.140

0.076

0.264

NSEI

p_pelagic

13

8

1

0.308

0.291

0.107

0.601

1.127

0.132

0.129

NSEI

p_pelagic

13

8

2

0.128

0.111

0.034

0.315

1.152

0.073

0.232

NSEI

p_pelagic

13

9

1

0.309

0.290

0.107

0.605

1.124

0.131

0.145

NSEI

p_pelagic

13

9

2

0.128

0.110

0.034

0.334

1.158

0.077

0.235

NSEI

p_pelagic

14

2

2

0.410

0.404

0.184

0.667

1.011

0.125

0.048

NSEO

p_black

14

3

1

0.918

0.945

0.688

0.993

1.017

0.083

-0.327

NSEO

p_pelagic

16

1

1

0.345

0.332

0.116

0.642

1.225

0.140

0.093

SSEO

p_pelagic

16

1

2

0.326

0.311

0.115

0.625

1.219

0.136

0.110

SSEO

p_pelagic

16

10

1

0.350

0.342

0.122

0.653

1.194

0.139

0.057

SSEO

p_pelagic

16

10

2

0.333

0.320

0.118

0.624

1.213

0.136

0.096

SSEO

p_pelagic

16

11

1

0.353

0.341

0.128

0.646

1.160

0.139

0.086

SSEO

p_pelagic

16

11

2

0.332

0.320

0.118

0.605

1.203

0.131

0.091

SSEO

p_pelagic

16

12

1

0.355

0.345

0.128

0.650

1.177

0.138

0.072

SSEO

p_pelagic

16

12

2

0.331

0.315

0.121

0.614

1.163

0.134

0.120

SSEO

p_pelagic

16

13

1

0.357

0.343

0.132

0.640

1.162

0.136

0.103

SSEO

p_pelagic

16

13

2

0.338

0.324

0.124

0.623

1.165

0.133

0.105

SSEO

p_pelagic

16

14

1

0.345

0.332

0.132

0.621

1.144

0.131

0.099

SSEO

p_pelagic

16

14

2

0.337

0.323

0.123

0.625

1.135

0.132

0.106

SSEO

p_pelagic

16

15

1

0.351

0.339

0.139

0.623

1.153

0.130

0.092

SSEO

p_pelagic

16

15

2

0.340

0.327

0.132

0.619

1.127

0.129

0.101

SSEO

p_pelagic

16

16

1

0.349

0.336

0.136

0.627

1.133

0.131

0.099

SSEO

p_pelagic

16

16

2

0.342

0.332

0.134

0.609

1.122

0.127

0.079

SSEO

p_pelagic

16

17

1

0.404

0.394

0.181

0.669

1.095

0.129

0.078

SSEO

p_pelagic

16

17

2

0.360

0.348

0.135

0.641

1.127

0.133

0.090

SSEO

p_pelagic

16

18

1

0.390

0.384

0.163

0.660

1.096

0.130

0.046

SSEO

p_pelagic

16

18

2

0.360

0.349

0.152

0.629

1.106

0.128

0.086

SSEO

p_pelagic

16

19

1

0.379

0.367

0.163

0.651

1.093

0.128

0.093

SSEO

p_pelagic

16

19

2

0.363

0.354

0.145

0.632

1.077

0.130

0.069

SSEO

p_pelagic

16

2

1

0.342

0.326

0.122

0.645

1.237

0.141

0.113

SSEO

p_pelagic

16

2

2

0.324

0.307

0.104

0.625

1.253

0.137

0.124

SSEO

p_pelagic

16

20

1

0.387

0.378

0.172

0.650

1.085

0.125

0.072

SSEO

p_pelagic

16

20

2

0.369

0.361

0.153

0.628

1.076

0.125

0.064

SSEO

p_pelagic

16

21

1

0.371

0.364

0.162

0.622

1.057

0.122

0.058

SSEO

p_pelagic

16

21

2

0.369

0.357

0.158

0.631

1.047

0.125

0.096

SSEO

p_pelagic

16

22

2

0.377

0.371

0.166

0.638

1.035

0.123

0.049

SSEO

p_pelagic

16

23

2

0.300

0.287

0.119

0.547

1.023

0.112

0.116

SSEO

p_pelagic

16

3

1

0.348

0.327

0.124

0.646

1.246

0.142

0.148

SSEO

p_pelagic

16

3

2

0.329

0.316

0.107

0.634

1.240

0.140

0.093

SSEO

p_dsr

16

30

1

0.419

0.419

0.393

0.447

1.015

0.014

0.000

SSEO

p_pelagic

16

4

1

0.344

0.329

0.122

0.632

1.245

0.140

0.107

SSEO

p_pelagic

16

4

2

0.325

0.311

0.105

0.626

1.202

0.136

0.103

SSEO

p_pelagic

16

5

1

0.347

0.333

0.123

0.639

1.214

0.140

0.100

SSEO

p_pelagic

16

5

2

0.325

0.308

0.107

0.611

1.190

0.136

0.125

SSEO

p_pelagic

16

6

1

0.346

0.330

0.118

0.634

1.208

0.141

0.114

SSEO

p_pelagic

16

6

2

0.326

0.317

0.100

0.615

1.185

0.136

0.066

SSEO

p_pelagic

16

7

1

0.345

0.332

0.119

0.631

1.204

0.138

0.094

SSEO

p_pelagic

16

7

2

0.331

0.314

0.117

0.619

1.223

0.135

0.126

SSEO

p_pelagic

16

8

1

0.347

0.332

0.124

0.639

1.214

0.138

0.109

SSEO

p_pelagic

16

8

2

0.330

0.315

0.110

0.626

1.195

0.137

0.110

SSEO

p_pelagic

16

9

1

0.352

0.341

0.120

0.649

1.187

0.140

0.079

SSEO

p_pelagic

16

9

2

0.329

0.313

0.116

0.618

1.193

0.133

0.121

SSEO

p_pelagic

3

1

1

0.546

0.558

0.262

0.780

1.011

0.135

-0.089

PWSI

p_black

8

40

1

0.523

0.422

0.341

0.969

1.143

0.220

0.459

afognak

p_black

8

48

1

0.537

0.493

0.411

0.969

1.043

0.148

0.298

afognak

Table Xb. True parameters with an Rhat in excess of 1.1

variable

area_n

year/species

mean

median

lower

upper

Rhat

sd

skew

year

area

beta2_pH

10

1

1.579

0.895

0.320

7.822

1.145

1.878

0.364

northeast

beta0_pH

11

4

0.860

0.897

-0.343

1.858

1.109

0.537

-0.069

CSEO

beta1_pH

11

1

3.488

3.441

2.436

4.653

1.123

0.568

0.083

CSEO

beta1_pH

11

4

1.926

1.918

0.786

3.104

1.103

0.579

0.014

CSEO

beta1_pelagic

11

1.185

1.158

0.711

1.824

1.142

0.263

0.103

CSEO

beta2_yellow

11

-6.842

-5.681

-20.016

-0.411

1.118

5.237

-0.222

CSEO

beta0_pH

13

4

0.478

0.588

-1.061

1.576

1.179

0.688

-0.160

NSEI

beta1_pH

13

4

1.494

1.401

0.346

3.045

1.178

0.688

0.135

NSEI

beta0_pelagic

13

-1.008

-1.103

-1.713

-0.004

1.341

0.475

0.200

NSEI

beta1_pelagic

13

2.810

2.891

1.477

3.843

1.325

0.614

-0.132

NSEI

beta2_pelagic

13

0.188

0.162

0.098

0.443

1.263

0.092

0.283

NSEI

beta3_pelagic

13

31.164

30.680

24.723

38.969

1.348

3.996

0.121

NSEI

beta0_pH

15

4

0.667

0.684

-0.291

1.576

1.561

0.459

-0.037

SSEI

beta1_pH

15

4

1.345

1.343

0.290

2.334

1.514

0.521

0.004

SSEI

beta1_pelagic

15

0.090

0.013

0.000

0.603

1.108

0.162

0.477

SSEI

beta0_pelagic

16

-0.719

-0.733

-1.524

0.026

1.695

0.455

0.031

SSEO

beta1_pelagic

16

1.858

1.878

0.827

3.034

1.663

0.629

-0.032

SSEO

beta2_pelagic

16

0.560

0.271

0.108

3.745

1.195

0.875

0.330

SSEO

beta3_pelagic

16

29.877

30.135

24.912

34.665

1.184

2.591

-0.100

SSEO

mu_beta2_pH

2

1

1.583

0.882

0.198

8.069

1.161

2.029

0.346

NG

beta0_pH

2

1

1.418

1.426

1.011

1.723

1.121

0.180

-0.044

NG

beta1_pH

2

1

1.545

1.535

1.037

2.134

1.104

0.282

0.036

NG

mu_beta0_pH

3

4

0.688

0.729

-0.392

1.582

1.141

0.512

-0.080

PWSI

tau_beta1_pH

3

4

5.104

3.020

0.247

21.889

1.144

7.293

0.286

PWSI

tau_beta2_pH

3

1

4.714

1.092

0.044

39.204

1.112

10.952

0.331

PWSI

beta0_pH

3

1

1.354

1.403

0.366

1.768

1.236

0.303

-0.162

PWSI

beta1_pH

3

1

1.865

1.785

1.108

3.437

1.207

0.514

0.156

PWSI

beta2_pH

6

1

1.361

0.523

0.126

8.858

1.146

2.224

0.377

SOKO2SAP

beta0_pH

8

1

-0.560

-0.539

-1.180

-0.033

1.175

0.294

-0.071

afognak

beta2_pH

8

1

1.024

0.402

0.150

6.871

1.256

1.714

0.363

afognak

beta3_pH

8

1

35.023

34.860

29.093

42.115

1.251

3.614

0.045

afognak

beta2_pH

9

1

1.365

0.610

0.158

8.275

1.112

1.990

0.379

eastside

lambda_H

15

0.019

0.014

0.003

0.057

1.249

0.047

0.107

SSEI

lambda_H

6

3.837

0.117

0.008

32.088

1.144

9.178

0.405

SOKO2SAP

Table Xc. Summary statistics describing the characteristics of the derived harvest and release estimates that are not converged relative to an Rhat of 1.01.

Statistic

non.BSAI.SOKO.Harv

BSAI.SOKO.Harv

non.BSAI.SOKO.Rel

BSAI.SOKO.Rel

n unconverged at Rhat = 1.01

636.00

292.00

1,186.00

358.00

Prop. of total

0.69

0.31

0.77

0.23

Prop. median less than 1

0.21

0.16

0.31

0.41

Prop. median less than 10

0.33

0.41

0.48

0.70

Prop. median less than 50

0.49

0.62

0.65

0.77

Median of the median values

53.28

21.26

11.16

1.88

Mean of the median values

584.32

627.73

332.50

261.96

Mean skew of posteriors

0.23

0.29

0.27

0.25

Table Xd. Summary statistics describing the characteristics of the derived harvest and release estimates that are not converged relative to an Rhat of 1.1.

Statistic

non.BSAI.SOKO.Harv

BSAI.SOKO.Harv

non.BSAI.SOKO.Rel

BSAI.SOKO.Rel

n unconverged at Rhat = 1.1

31.00

3.00

68.00

35.00

Prop. of total

0.91

0.09

0.66

0.34

Prop. median less than 1

0.13

1.00

0.71

0.63

Prop. median less than 10

0.16

1.00

0.76

0.91

Prop. median less than 50

0.19

1.00

0.82

0.91

Median of the median values

295.80

0.09

0.01

0.28

Mean of the median values

443.61

0.36

51.42

14.85

Mean skew of posteriors

0.25

0.21

0.16

0.16

Table Xe. All parameters and derived quantities that are not converged relative to an Rhat of 1.1.

variable

Freq

beta0_pelagic

2

beta0_pH

5

beta1_pelagic

3

beta1_pH

4

beta2_pelagic

2

beta2_pH

4

beta2_yellow

1

beta3_pelagic

2

beta3_pH

1

Bp_ay

5

Bp_ayu

4

eps_bl

1

eps_pel

6

eps_pH

11

H_ay

2

H_ayg

1

H_ayu

1

Hb_ay

6

Hb_ayg

2

Hb_ayu

3

Ho_ay

2

Ho_ayg

2

Ho_ayu

3

Hp_ay

12

Hp_ayg

3

Hp_ayu

7

Hs_ayu

1

Htrend_ay

4

Hy_ayu

1

lambda_H

2

mu_beta0_pH

1

mu_beta2_pH

1

p_black

1

p_pelagic

62

p_slope

1

pH

5

R_ay

2

R_ayg

6

R_ayu

8

Rb_ay

1

Rb_ayg

5

Rb_ayu

1

Rd_ay

1

Rd_ayg

2

Rd_ayu

1

re_black

1

re_pelagic

19

re_pH

16

Ro_ay

18

Ro_ayg

13

Ro_ayu

33

Rp_ay

2

Rp_ayg

6

Rp_ayu

2

Rs_ay

3

Rs_ayg

9

Rs_ayu

16

Ry_ay

11

Ry_ayg

4

Ry_ayu

23

tau_beta1_pH

1

tau_beta2_pH

1

Parameter estimates:

Summary Table of Parameter Estimates
Parameter mean sd Lower_CI Median Upper_CI
mu_bc_H[1] -0.020 0.067 -0.145 -0.023 0.124
mu_bc_H[2] 0.078 0.054 -0.018 0.076 0.190
mu_bc_H[3] -0.267 0.068 -0.395 -0.266 -0.134
mu_bc_H[4] -0.921 0.179 -1.280 -0.921 -0.576
mu_bc_H[5] 1.308 1.064 0.140 1.071 3.942
mu_bc_H[6] -2.376 0.313 -2.999 -2.372 -1.773
mu_bc_H[7] -0.456 0.103 -0.665 -0.455 -0.262
mu_bc_H[8] 0.253 0.356 -0.326 0.222 0.984
mu_bc_H[9] -0.244 0.126 -0.481 -0.248 0.015
mu_bc_H[10] 0.058 0.067 -0.069 0.056 0.198
mu_bc_H[11] -0.124 0.037 -0.197 -0.124 -0.050
mu_bc_H[12] -0.254 0.105 -0.476 -0.251 -0.053
mu_bc_H[13] -0.157 0.072 -0.290 -0.159 -0.020
mu_bc_H[14] -0.311 0.087 -0.491 -0.309 -0.147
mu_bc_H[15] -0.338 0.048 -0.433 -0.336 -0.242
mu_bc_H[16] -0.161 0.371 -0.805 -0.186 0.649
mu_bc_R[1] 1.041 0.125 0.799 1.040 1.279
mu_bc_R[2] 1.014 0.076 0.854 1.017 1.156
mu_bc_R[3] 0.226 0.092 0.044 0.226 0.404
mu_bc_R[4] 0.676 0.197 0.271 0.683 1.049
mu_bc_R[5] 1.686 0.367 0.970 1.680 2.433
mu_bc_R[6] -1.404 0.349 -2.096 -1.413 -0.722
mu_bc_R[7] 0.505 0.169 0.158 0.509 0.826
mu_bc_R[8] 0.657 0.149 0.378 0.657 0.954
mu_bc_R[9] 0.403 0.163 0.062 0.409 0.708
mu_bc_R[10] 0.899 0.105 0.693 0.901 1.105
mu_bc_R[11] 1.106 0.074 0.962 1.105 1.257
mu_bc_R[12] 0.874 0.165 0.564 0.872 1.196
mu_bc_R[13] 1.037 0.091 0.860 1.038 1.212
mu_bc_R[14] 0.855 0.119 0.624 0.857 1.090
mu_bc_R[15] 0.798 0.095 0.606 0.799 0.988
mu_bc_R[16] 1.357 0.115 1.126 1.361 1.582
tau_pH[1] 2.432 0.182 2.090 2.423 2.794
tau_pH[2] 1.611 0.150 1.337 1.610 1.914
tau_pH[3] 1.194 0.110 0.993 1.191 1.416
tau_pH[4] 1.015 0.138 0.770 1.007 1.307
tau_pH[5] 0.427 0.076 0.297 0.420 0.592
beta0_pH[1,1] 0.581 0.171 0.219 0.587 0.898
beta0_pH[2,1] 1.418 0.180 1.011 1.426 1.723
beta0_pH[3,1] 1.354 0.303 0.366 1.403 1.768
beta0_pH[4,1] 1.572 0.251 0.997 1.593 2.004
beta0_pH[5,1] -0.735 0.376 -1.589 -0.710 -0.092
beta0_pH[6,1] -0.476 0.452 -1.732 -0.354 0.090
beta0_pH[7,1] 0.757 0.011 0.736 0.757 0.778
beta0_pH[8,1] -0.560 0.294 -1.180 -0.539 -0.033
beta0_pH[9,1] -0.249 0.279 -0.868 -0.212 0.190
beta0_pH[10,1] 0.278 0.190 -0.094 0.287 0.653
beta0_pH[11,1] -0.032 0.337 -0.662 -0.062 0.719
beta0_pH[12,1] 0.632 0.198 0.216 0.642 1.000
beta0_pH[13,1] -0.654 0.150 -0.955 -0.651 -0.376
beta0_pH[14,1] 0.004 0.176 -0.363 0.006 0.339
beta0_pH[15,1] 0.041 0.185 -0.348 0.051 0.374
beta0_pH[16,1] -0.914 0.220 -1.346 -0.915 -0.469
beta0_pH[1,2] 2.563 0.084 2.404 2.560 2.734
beta0_pH[2,2] 3.079 0.022 3.036 3.080 3.123
beta0_pH[3,2] 3.626 0.042 3.544 3.625 3.711
beta0_pH[4,2] 2.937 0.019 2.901 2.936 2.974
beta0_pH[5,2] 4.528 1.488 2.713 4.229 8.362
beta0_pH[6,2] 3.320 0.163 3.005 3.316 3.642
beta0_pH[7,2] 2.033 0.049 1.941 2.032 2.133
beta0_pH[8,2] 2.820 0.065 2.697 2.819 2.953
beta0_pH[9,2] 3.954 0.210 3.569 3.947 4.382
beta0_pH[10,2] 3.996 0.130 3.746 3.993 4.258
beta0_pH[11,2] 1.885 0.009 1.869 1.885 1.903
beta0_pH[12,2] 1.796 0.069 1.662 1.794 1.936
beta0_pH[13,2] 2.496 0.026 2.446 2.496 2.548
beta0_pH[14,2] 1.803 0.024 1.759 1.803 1.850
beta0_pH[15,2] 2.545 0.017 2.512 2.545 2.578
beta0_pH[16,2] 2.765 0.024 2.718 2.765 2.812
beta0_pH[1,3] 1.075 0.401 0.151 1.146 1.679
beta0_pH[2,3] 2.126 0.018 2.092 2.126 2.160
beta0_pH[3,3] 2.692 0.024 2.645 2.692 2.740
beta0_pH[4,3] 2.980 0.034 2.915 2.980 3.046
beta0_pH[5,3] 3.495 2.442 0.911 2.854 9.990
beta0_pH[6,3] 1.382 0.381 0.721 1.362 2.187
beta0_pH[7,3] 0.859 0.051 0.760 0.861 0.958
beta0_pH[8,3] 0.249 0.301 -0.469 0.291 0.653
beta0_pH[9,3] -2.299 1.062 -4.732 -2.185 -0.501
beta0_pH[10,3] -0.427 0.835 -2.414 -0.327 0.891
beta0_pH[11,3] 1.157 0.486 0.132 1.175 2.049
beta0_pH[12,3] 1.009 0.651 -0.366 0.999 2.261
beta0_pH[13,3] 1.047 0.481 -0.004 1.072 1.909
beta0_pH[14,3] 0.620 0.413 -0.277 0.651 1.364
beta0_pH[15,3] 0.774 0.414 -0.087 0.790 1.528
beta0_pH[16,3] 1.185 0.634 -0.095 1.181 2.465
beta0_pH[11,4] 0.860 0.537 -0.343 0.897 1.858
beta0_pH[12,4] 0.883 0.716 -0.662 0.885 2.246
beta0_pH[13,4] 0.478 0.688 -1.061 0.588 1.576
beta0_pH[14,4] 0.392 0.466 -0.630 0.429 1.240
beta0_pH[15,4] 0.667 0.459 -0.291 0.684 1.576
beta0_pH[16,4] 0.855 0.632 -0.495 0.889 2.047
beta0_pH[11,5] 2.971 0.145 2.698 2.968 3.264
beta0_pH[12,5] -0.920 0.518 -2.018 -0.897 0.014
beta0_pH[13,5] 4.902 0.422 4.144 4.885 5.794
beta0_pH[14,5] 3.894 0.364 3.233 3.881 4.669
beta0_pH[15,5] 2.285 0.124 2.049 2.282 2.540
beta0_pH[16,5] 0.253 0.065 0.161 0.240 0.410
beta1_pH[1,1] 2.476 0.464 1.704 2.420 3.557
beta1_pH[2,1] 1.545 0.282 1.037 1.535 2.134
beta1_pH[3,1] 1.865 0.514 1.108 1.785 3.437
beta1_pH[4,1] 1.893 0.417 1.214 1.838 2.882
beta1_pH[5,1] 2.242 0.484 1.407 2.202 3.346
beta1_pH[6,1] 2.027 0.757 1.065 1.823 4.104
beta1_pH[7,1] 0.000 0.000 0.000 0.000 0.000
beta1_pH[8,1] 2.228 0.777 1.136 2.084 4.126
beta1_pH[9,1] 1.968 0.600 1.188 1.847 3.688
beta1_pH[10,1] 2.098 0.362 1.456 2.084 2.878
beta1_pH[11,1] 3.488 0.568 2.436 3.441 4.653
beta1_pH[12,1] 2.894 0.292 2.341 2.885 3.469
beta1_pH[13,1] 3.451 0.239 3.018 3.438 3.940
beta1_pH[14,1] 2.857 0.330 2.263 2.840 3.570
beta1_pH[15,1] 2.414 0.252 1.969 2.405 2.941
beta1_pH[16,1] 4.087 0.303 3.472 4.093 4.689
beta1_pH[1,2] 0.000 0.000 0.000 0.000 0.000
beta1_pH[2,2] 0.000 0.000 0.000 0.000 0.000
beta1_pH[3,2] 0.000 0.000 0.000 0.000 0.000
beta1_pH[4,2] 0.000 0.000 0.000 0.000 0.000
beta1_pH[5,2] 0.000 0.000 0.000 0.000 0.000
beta1_pH[6,2] 0.000 0.000 0.000 0.000 0.000
beta1_pH[7,2] 0.000 0.000 0.000 0.000 0.000
beta1_pH[8,2] 0.000 0.000 0.000 0.000 0.000
beta1_pH[9,2] 0.000 0.000 0.000 0.000 0.000
beta1_pH[10,2] 0.000 0.000 0.000 0.000 0.000
beta1_pH[11,2] 0.000 0.000 0.000 0.000 0.000
beta1_pH[12,2] 0.000 0.000 0.000 0.000 0.000
beta1_pH[13,2] 0.000 0.000 0.000 0.000 0.000
beta1_pH[14,2] 0.000 0.000 0.000 0.000 0.000
beta1_pH[15,2] 0.000 0.000 0.000 0.000 0.000
beta1_pH[16,2] 0.000 0.000 0.000 0.000 0.000
beta1_pH[1,3] 1.969 1.174 0.049 1.835 5.116
beta1_pH[2,3] 0.000 0.000 0.000 0.000 0.000
beta1_pH[3,3] 0.000 0.000 0.000 0.000 0.000
beta1_pH[4,3] 0.000 0.000 0.000 0.000 0.000
beta1_pH[5,3] 0.000 0.000 0.000 0.000 0.000
beta1_pH[6,3] 0.000 0.000 0.000 0.000 0.000
beta1_pH[7,3] 0.000 0.000 0.000 0.000 0.000
beta1_pH[8,3] 4.735 2.537 0.619 4.224 11.082
beta1_pH[9,3] 7.669 2.717 3.360 7.305 13.825
beta1_pH[10,3] 6.277 2.571 2.717 5.789 12.506
beta1_pH[11,3] 2.127 0.523 1.110 2.119 3.177
beta1_pH[12,3] 3.214 0.669 1.933 3.203 4.654
beta1_pH[13,3] 1.929 0.517 0.921 1.912 3.025
beta1_pH[14,3] 2.148 0.469 1.273 2.136 3.140
beta1_pH[15,3] 2.039 0.484 1.123 2.026 3.018
beta1_pH[16,3] 3.695 0.685 2.329 3.706 5.110
beta1_pH[11,4] 1.926 0.579 0.786 1.918 3.104
beta1_pH[12,4] 3.189 0.744 1.798 3.156 4.763
beta1_pH[13,4] 1.494 0.688 0.346 1.401 3.045
beta1_pH[14,4] 1.976 0.538 0.960 1.954 3.137
beta1_pH[15,4] 1.345 0.521 0.290 1.343 2.334
beta1_pH[16,4] 3.536 0.695 2.186 3.500 5.032
beta1_pH[11,5] 0.000 0.000 0.000 0.000 0.000
beta1_pH[12,5] 0.000 0.000 0.000 0.000 0.000
beta1_pH[13,5] 0.000 0.000 0.000 0.000 0.000
beta1_pH[14,5] 0.000 0.000 0.000 0.000 0.000
beta1_pH[15,5] 0.000 0.000 0.000 0.000 0.000
beta1_pH[16,5] 0.000 0.000 0.000 0.000 0.000
beta2_pH[1,1] 1.379 0.890 0.344 1.085 3.377
beta2_pH[2,1] 1.767 0.929 0.386 1.619 3.430
beta2_pH[3,1] 1.383 0.863 0.159 1.274 2.946
beta2_pH[4,1] 1.171 0.797 0.221 0.932 2.895
beta2_pH[5,1] 2.298 2.452 0.293 1.549 9.821
beta2_pH[6,1] 1.361 2.224 0.126 0.523 8.858
beta2_pH[7,1] 0.000 0.000 0.000 0.000 0.000
beta2_pH[8,1] 1.024 1.714 0.150 0.402 6.871
beta2_pH[9,1] 1.365 1.990 0.158 0.610 8.275
beta2_pH[10,1] 1.579 1.878 0.320 0.895 7.822
beta2_pH[11,1] 0.379 0.201 0.172 0.325 0.967
beta2_pH[12,1] 1.573 0.807 0.614 1.336 3.672
beta2_pH[13,1] 1.032 0.446 0.545 0.918 2.306
beta2_pH[14,1] 1.093 0.626 0.387 0.933 3.134
beta2_pH[15,1] 0.991 0.559 0.380 0.852 2.653
beta2_pH[16,1] 1.003 0.310 0.613 0.932 1.815
beta2_pH[1,2] 0.000 0.000 0.000 0.000 0.000
beta2_pH[2,2] 0.000 0.000 0.000 0.000 0.000
beta2_pH[3,2] 0.000 0.000 0.000 0.000 0.000
beta2_pH[4,2] 0.000 0.000 0.000 0.000 0.000
beta2_pH[5,2] 0.000 0.000 0.000 0.000 0.000
beta2_pH[6,2] 0.000 0.000 0.000 0.000 0.000
beta2_pH[7,2] 0.000 0.000 0.000 0.000 0.000
beta2_pH[8,2] 0.000 0.000 0.000 0.000 0.000
beta2_pH[9,2] 0.000 0.000 0.000 0.000 0.000
beta2_pH[10,2] 0.000 0.000 0.000 0.000 0.000
beta2_pH[11,2] 0.000 0.000 0.000 0.000 0.000
beta2_pH[12,2] 0.000 0.000 0.000 0.000 0.000
beta2_pH[13,2] 0.000 0.000 0.000 0.000 0.000
beta2_pH[14,2] 0.000 0.000 0.000 0.000 0.000
beta2_pH[15,2] 0.000 0.000 0.000 0.000 0.000
beta2_pH[16,2] 0.000 0.000 0.000 0.000 0.000
beta2_pH[1,3] 6.654 6.578 0.154 4.905 22.856
beta2_pH[2,3] 0.000 0.000 0.000 0.000 0.000
beta2_pH[3,3] 0.000 0.000 0.000 0.000 0.000
beta2_pH[4,3] 0.000 0.000 0.000 0.000 0.000
beta2_pH[5,3] 0.000 0.000 0.000 0.000 0.000
beta2_pH[6,3] 0.000 0.000 0.000 0.000 0.000
beta2_pH[7,3] 0.000 0.000 0.000 0.000 0.000
beta2_pH[8,3] 1.539 1.531 0.118 1.055 6.375
beta2_pH[9,3] 0.348 0.883 0.088 0.195 1.404
beta2_pH[10,3] 0.481 0.917 0.109 0.278 2.172
beta2_pH[11,3] -6.286 4.856 -18.752 -4.947 -0.724
beta2_pH[12,3] -6.934 4.737 -19.149 -5.890 -1.220
beta2_pH[13,3] -6.427 4.722 -18.635 -5.015 -1.076
beta2_pH[14,3] -6.393 4.628 -18.457 -4.979 -1.285
beta2_pH[15,3] -6.533 4.667 -18.559 -5.273 -1.040
beta2_pH[16,3] -6.522 4.577 -18.808 -5.156 -1.507
beta2_pH[11,4] -5.273 3.948 -14.736 -4.077 -0.752
beta2_pH[12,4] -5.660 4.015 -15.711 -4.542 -1.027
beta2_pH[13,4] -5.380 3.967 -15.173 -4.143 -0.814
beta2_pH[14,4] -5.521 3.910 -15.256 -4.414 -1.069
beta2_pH[15,4] -5.570 3.974 -15.302 -4.470 -0.844
beta2_pH[16,4] -5.610 3.748 -15.015 -4.512 -1.377
beta2_pH[11,5] 0.000 0.000 0.000 0.000 0.000
beta2_pH[12,5] 0.000 0.000 0.000 0.000 0.000
beta2_pH[13,5] 0.000 0.000 0.000 0.000 0.000
beta2_pH[14,5] 0.000 0.000 0.000 0.000 0.000
beta2_pH[15,5] 0.000 0.000 0.000 0.000 0.000
beta2_pH[16,5] 0.000 0.000 0.000 0.000 0.000
beta3_pH[1,1] 34.866 1.083 32.828 34.886 37.279
beta3_pH[2,1] 32.368 1.303 29.786 32.368 34.953
beta3_pH[3,1] 32.919 1.469 29.927 32.917 35.868
beta3_pH[4,1] 32.642 1.829 29.677 32.513 36.433
beta3_pH[5,1] 27.694 1.235 26.131 27.515 31.175
beta3_pH[6,1] 33.121 3.115 26.663 32.974 40.789
beta3_pH[7,1] 25.000 0.000 25.000 25.000 25.000
beta3_pH[8,1] 35.023 3.614 29.093 34.860 42.115
beta3_pH[9,1] 32.569 2.569 28.179 32.199 38.639
beta3_pH[10,1] 32.358 1.239 30.082 32.315 34.976
beta3_pH[11,1] 31.163 1.072 29.273 31.091 33.275
beta3_pH[12,1] 30.790 0.467 29.810 30.803 31.680
beta3_pH[13,1] 31.423 0.539 30.394 31.417 32.485
beta3_pH[14,1] 32.494 0.742 31.199 32.455 34.039
beta3_pH[15,1] 30.915 0.800 29.342 30.918 32.537
beta3_pH[16,1] 30.106 0.354 29.438 30.099 30.817
beta3_pH[1,2] 34.543 5.416 25.577 34.480 43.495
beta3_pH[2,2] 34.703 5.423 25.474 34.852 43.517
beta3_pH[3,2] 30.405 7.179 18.676 30.252 42.362
beta3_pH[4,2] 30.306 7.166 18.599 30.263 42.318
beta3_pH[5,2] 30.506 7.297 18.568 30.500 42.388
beta3_pH[6,2] 30.741 7.220 18.558 30.755 42.464
beta3_pH[7,2] 25.000 0.000 25.000 25.000 25.000
beta3_pH[8,2] 25.000 0.000 25.000 25.000 25.000
beta3_pH[9,2] 25.000 0.000 25.000 25.000 25.000
beta3_pH[10,2] 25.000 0.000 25.000 25.000 25.000
beta3_pH[11,2] 24.000 0.000 24.000 24.000 24.000
beta3_pH[12,2] 24.000 0.000 24.000 24.000 24.000
beta3_pH[13,2] 24.000 0.000 24.000 24.000 24.000
beta3_pH[14,2] 24.000 0.000 24.000 24.000 24.000
beta3_pH[15,2] 24.000 0.000 24.000 24.000 24.000
beta3_pH[16,2] 24.000 0.000 24.000 24.000 24.000
beta3_pH[1,3] 38.305 3.199 30.847 39.597 42.481
beta3_pH[2,3] 34.469 5.455 25.475 34.489 43.516
beta3_pH[3,3] 30.377 7.135 18.637 30.394 42.248
beta3_pH[4,3] 30.500 7.138 18.784 30.477 42.510
beta3_pH[5,3] 25.000 0.000 25.000 25.000 25.000
beta3_pH[6,3] 25.000 0.000 25.000 25.000 25.000
beta3_pH[7,3] 25.000 0.000 25.000 25.000 25.000
beta3_pH[8,3] 42.028 1.373 38.275 42.082 43.903
beta3_pH[9,3] 36.204 1.900 34.528 35.442 41.810
beta3_pH[10,3] 37.755 2.448 34.672 37.071 43.421
beta3_pH[11,3] 40.696 1.069 38.800 40.653 42.736
beta3_pH[12,3] 41.992 0.742 40.319 42.244 42.859
beta3_pH[13,3] 41.399 0.730 40.009 41.409 42.829
beta3_pH[14,3] 40.477 0.675 39.128 40.494 41.730
beta3_pH[15,3] 40.629 0.650 39.194 40.641 41.893
beta3_pH[16,3] 41.846 0.391 41.044 41.907 42.622
beta3_pH[11,4] 40.616 1.187 38.359 40.494 42.981
beta3_pH[12,4] 41.857 0.787 40.223 42.111 42.950
beta3_pH[13,4] 42.268 1.190 39.837 42.566 43.895
beta3_pH[14,4] 40.479 0.734 38.998 40.524 41.802
beta3_pH[15,4] 41.100 1.147 38.360 41.056 43.418
beta3_pH[16,4] 41.743 0.423 40.928 41.785 42.567
beta3_pH[11,5] 25.000 0.000 25.000 25.000 25.000
beta3_pH[12,5] 25.000 0.000 25.000 25.000 25.000
beta3_pH[13,5] 25.000 0.000 25.000 25.000 25.000
beta3_pH[14,5] 25.000 0.000 25.000 25.000 25.000
beta3_pH[15,5] 25.000 0.000 25.000 25.000 25.000
beta3_pH[16,5] 25.000 0.000 25.000 25.000 25.000
beta4_pH[1,1] 0.264 0.021 0.226 0.264 0.307
beta4_pH[2,1] 0.211 0.005 0.201 0.211 0.221
beta4_pH[3,1] 0.173 0.007 0.159 0.173 0.186
beta4_pH[4,1] 0.294 0.012 0.271 0.294 0.316
beta4_pH[5,1] 0.419 0.277 0.086 0.385 0.922
beta4_pH[6,1] 0.784 0.402 0.325 0.733 1.445
beta4_pH[7,1] 0.858 0.340 0.355 0.809 1.646
beta4_pH[8,1] 0.159 0.069 0.046 0.153 0.309
beta4_pH[9,1] 0.211 0.084 0.062 0.208 0.388
beta4_pH[10,1] 0.001 0.001 0.000 0.001 0.003
beta4_pH[11,1] 0.447 0.012 0.424 0.447 0.470
beta4_pH[12,1] 0.641 0.022 0.598 0.641 0.686
beta4_pH[13,1] 0.631 0.013 0.605 0.631 0.657
beta4_pH[14,1] 0.670 0.030 0.612 0.669 0.731
beta4_pH[15,1] 0.796 0.021 0.756 0.797 0.838
beta4_pH[16,1] 1.218 0.042 1.141 1.216 1.306
beta4_pH[1,2] 0.308 0.094 0.136 0.300 0.507
beta4_pH[2,2] 0.514 0.023 0.470 0.513 0.561
beta4_pH[3,2] 0.446 0.012 0.422 0.446 0.469
beta4_pH[4,2] 0.654 0.038 0.581 0.653 0.730
beta4_pH[5,2] 0.733 0.918 0.112 0.530 2.733
beta4_pH[6,2] 0.916 0.795 0.202 0.737 2.876
beta4_pH[7,2] 0.894 0.446 0.357 0.810 1.914
beta4_pH[8,2] 0.849 0.199 0.494 0.837 1.265
beta4_pH[9,2] 1.195 1.235 0.306 0.886 4.042
beta4_pH[10,2] 0.457 0.051 0.361 0.457 0.559
beta4_pH[11,2] 1.270 0.049 1.177 1.270 1.368
beta4_pH[12,2] 0.581 0.114 0.373 0.573 0.829
beta4_pH[13,2] 0.945 0.045 0.859 0.944 1.037
beta4_pH[14,2] 0.931 0.064 0.803 0.931 1.057
beta4_pH[15,2] 1.011 0.025 0.962 1.011 1.061
beta4_pH[16,2] 1.407 0.077 1.259 1.405 1.559
beta4_pH[1,3] 0.288 0.086 0.135 0.283 0.464
beta4_pH[2,3] 0.052 0.017 0.020 0.051 0.087
beta4_pH[3,3] 0.053 0.010 0.034 0.052 0.072
beta4_pH[4,3] 0.173 0.031 0.113 0.172 0.234
beta4_pH[5,3] 0.777 0.954 0.122 0.563 2.769
beta4_pH[6,3] 0.854 0.679 0.199 0.702 2.493
beta4_pH[7,3] 0.879 0.398 0.346 0.806 1.867
beta4_pH[8,3] 0.418 0.194 0.110 0.396 0.852
beta4_pH[9,3] 0.939 0.937 0.177 0.708 2.954
beta4_pH[10,3] 0.061 0.041 0.003 0.056 0.155
beta4_pH[11,3] 0.364 0.016 0.333 0.363 0.396
beta4_pH[12,3] 0.440 0.041 0.361 0.438 0.524
beta4_pH[13,3] 0.839 0.037 0.771 0.838 0.914
beta4_pH[14,3] 0.825 0.066 0.703 0.821 0.959
beta4_pH[15,3] 0.636 0.017 0.604 0.635 0.668
beta4_pH[16,3] 0.638 0.035 0.573 0.637 0.709
beta4_pH[1,4] 0.624 0.330 0.212 0.559 1.445
beta4_pH[2,4] 0.808 0.875 0.123 0.585 2.852
beta4_pH[3,4] 0.850 0.862 0.113 0.607 3.084
beta4_pH[4,4] 0.725 0.449 0.201 0.617 1.829
beta4_pH[5,4] 0.762 0.881 0.119 0.559 2.531
beta4_pH[6,4] 0.941 1.037 0.204 0.745 2.859
beta4_pH[7,4] 0.896 0.439 0.349 0.825 1.897
beta4_pH[8,4] 0.750 0.632 0.165 0.590 2.315
beta4_pH[9,4] 0.763 0.750 0.135 0.569 2.544
beta4_pH[10,4] 1.454 2.910 0.008 0.404 10.628
beta4_pH[11,4] 0.444 0.017 0.412 0.444 0.477
beta4_pH[12,4] 0.562 0.053 0.458 0.561 0.667
beta4_pH[13,4] 0.925 0.036 0.856 0.924 0.996
beta4_pH[14,4] 0.919 0.073 0.778 0.918 1.070
beta4_pH[15,4] 0.713 0.020 0.674 0.713 0.754
beta4_pH[16,4] 0.686 0.042 0.608 0.684 0.773
beta4_pH[1,5] 0.630 0.333 0.214 0.560 1.383
beta4_pH[2,5] 0.813 0.890 0.108 0.592 2.822
beta4_pH[3,5] 0.855 0.978 0.117 0.588 3.166
beta4_pH[4,5] 0.735 0.490 0.209 0.629 1.967
beta4_pH[5,5] 0.761 0.784 0.105 0.559 2.747
beta4_pH[6,5] 0.919 0.991 0.214 0.727 2.653
beta4_pH[7,5] 0.896 0.439 0.340 0.812 1.900
beta4_pH[8,5] 0.746 0.650 0.164 0.603 2.140
beta4_pH[9,5] 0.760 0.765 0.126 0.560 2.687
beta4_pH[10,5] 1.446 2.986 0.007 0.360 11.347
beta4_pH[11,5] 0.602 0.053 0.507 0.600 0.711
beta4_pH[12,5] 0.583 0.164 0.307 0.564 0.951
beta4_pH[13,5] 0.773 0.063 0.661 0.769 0.906
beta4_pH[14,5] 0.842 0.118 0.633 0.835 1.095
beta4_pH[15,5] 1.145 0.072 1.012 1.144 1.290
beta4_pH[16,5] 18.618 4.019 10.580 18.909 24.686
beta0_pelagic[1] 2.269 0.037 2.201 2.268 2.343
beta0_pelagic[2] 1.537 0.007 1.524 1.537 1.551
beta0_pelagic[3] 0.177 0.275 -0.543 0.249 0.519
beta0_pelagic[4] 0.231 0.197 -0.172 0.243 0.574
beta0_pelagic[5] 1.446 0.156 1.156 1.441 1.777
beta0_pelagic[6] 1.354 0.061 1.239 1.353 1.483
beta0_pelagic[7] 1.480 0.039 1.408 1.478 1.558
beta0_pelagic[8] 1.673 0.043 1.591 1.672 1.761
beta0_pelagic[9] 2.794 0.066 2.670 2.793 2.934
beta0_pelagic[10] 2.603 0.031 2.544 2.603 2.664
beta0_pelagic[11] -0.164 0.181 -0.573 -0.153 0.142
beta0_pelagic[12] 1.801 0.036 1.731 1.800 1.873
beta0_pelagic[13] -1.008 0.475 -1.713 -1.103 -0.004
beta0_pelagic[14] -0.273 0.183 -0.714 -0.255 0.033
beta0_pelagic[15] -0.210 0.034 -0.301 -0.196 -0.181
beta0_pelagic[16] -0.719 0.455 -1.524 -0.733 0.026
beta1_pelagic[1] 0.000 0.000 0.000 0.000 0.000
beta1_pelagic[2] 0.000 0.000 0.000 0.000 0.000
beta1_pelagic[3] 0.785 0.314 0.375 0.716 1.708
beta1_pelagic[4] 0.814 0.218 0.429 0.811 1.263
beta1_pelagic[5] 0.000 0.000 0.000 0.000 0.000
beta1_pelagic[6] 0.000 0.000 0.000 0.000 0.000
beta1_pelagic[7] 0.000 0.000 0.000 0.000 0.000
beta1_pelagic[8] 0.000 0.000 0.000 0.000 0.000
beta1_pelagic[9] 0.000 0.000 0.000 0.000 0.000
beta1_pelagic[10] 0.000 0.000 0.000 0.000 0.000
beta1_pelagic[11] 1.185 0.263 0.711 1.158 1.824
beta1_pelagic[12] 0.000 0.000 0.000 0.000 0.000
beta1_pelagic[13] 2.810 0.614 1.477 2.891 3.843
beta1_pelagic[14] 1.813 0.474 1.100 1.713 2.969
beta1_pelagic[15] 0.090 0.162 0.000 0.013 0.603
beta1_pelagic[16] 1.858 0.629 0.827 1.878 3.034
beta2_pelagic[1] 0.000 0.000 0.000 0.000 0.000
beta2_pelagic[2] 0.000 0.000 0.000 0.000 0.000
beta2_pelagic[3] 6.006 5.950 0.171 4.272 21.267
beta2_pelagic[4] 7.477 6.054 0.408 6.083 22.113
beta2_pelagic[5] 0.000 0.000 0.000 0.000 0.000
beta2_pelagic[6] 0.000 0.000 0.000 0.000 0.000
beta2_pelagic[7] 0.000 0.000 0.000 0.000 0.000
beta2_pelagic[8] 0.000 0.000 0.000 0.000 0.000
beta2_pelagic[9] 0.000 0.000 0.000 0.000 0.000
beta2_pelagic[10] 0.000 0.000 0.000 0.000 0.000
beta2_pelagic[11] 1.661 1.488 0.220 1.166 5.959
beta2_pelagic[12] 0.000 0.000 0.000 0.000 0.000
beta2_pelagic[13] 0.188 0.092 0.098 0.162 0.443
beta2_pelagic[14] 0.813 0.828 0.199 0.547 3.195
beta2_pelagic[15] 2.639 1.902 0.138 2.257 7.300
beta2_pelagic[16] 0.560 0.875 0.108 0.271 3.745
beta3_pelagic[1] 28.687 3.868 22.287 28.689 35.101
beta3_pelagic[2] 28.753 3.890 22.338 28.572 35.155
beta3_pelagic[3] 27.911 2.682 23.055 28.252 32.184
beta3_pelagic[4] 26.074 1.241 23.266 26.284 28.643
beta3_pelagic[5] 29.312 6.688 18.659 28.966 40.365
beta3_pelagic[6] 29.541 6.650 18.545 29.556 40.475
beta3_pelagic[7] 29.533 6.672 18.537 29.438 40.461
beta3_pelagic[8] 29.454 6.593 18.557 29.231 40.364
beta3_pelagic[9] 29.598 6.659 18.584 29.574 40.441
beta3_pelagic[10] 29.386 6.480 18.580 29.337 40.245
beta3_pelagic[11] 30.042 1.325 27.596 29.974 32.875
beta3_pelagic[12] 30.000 0.000 30.000 30.000 30.000
beta3_pelagic[13] 31.164 3.996 24.723 30.680 38.969
beta3_pelagic[14] 34.470 1.743 31.516 34.239 38.336
beta3_pelagic[15] 33.391 5.906 23.536 33.778 41.723
beta3_pelagic[16] 29.877 2.591 24.912 30.135 34.665
beta0_yellow[1] -0.483 0.080 -0.638 -0.483 -0.324
beta0_yellow[2] 0.544 0.079 0.368 0.550 0.682
beta0_yellow[3] -0.220 0.133 -0.495 -0.221 0.013
beta0_yellow[4] 0.877 0.252 0.236 0.939 1.185
beta0_yellow[5] 0.248 0.335 -0.458 0.272 0.834
beta0_yellow[6] 0.822 0.107 0.614 0.825 1.025
beta0_yellow[7] 0.993 0.095 0.809 0.993 1.176
beta0_yellow[8] 1.023 0.092 0.846 1.022 1.205
beta0_yellow[9] 0.460 0.100 0.262 0.459 0.657
beta0_yellow[10] 0.540 0.072 0.399 0.538 0.683
beta0_yellow[11] -0.194 0.122 -0.460 -0.184 0.020
beta0_yellow[12] -1.081 0.094 -1.267 -1.079 -0.904
beta0_yellow[13] -1.362 0.185 -1.788 -1.345 -1.064
beta0_yellow[14] -0.248 0.142 -0.546 -0.234 -0.011
beta0_yellow[15] -1.254 0.139 -1.604 -1.236 -1.037
beta0_yellow[16] -0.244 0.093 -0.428 -0.244 -0.066
beta1_yellow[1] 0.000 0.000 0.000 0.000 0.000
beta1_yellow[2] 1.011 0.231 0.584 1.003 1.474
beta1_yellow[3] 0.483 0.246 0.020 0.497 0.936
beta1_yellow[4] 1.344 0.619 0.637 1.182 3.268
beta1_yellow[5] 0.000 0.000 0.000 0.000 0.000
beta1_yellow[6] 0.000 0.000 0.000 0.000 0.000
beta1_yellow[7] 0.000 0.000 0.000 0.000 0.000
beta1_yellow[8] 0.000 0.000 0.000 0.000 0.000
beta1_yellow[9] 0.000 0.000 0.000 0.000 0.000
beta1_yellow[10] 0.000 0.000 0.000 0.000 0.000
beta1_yellow[11] 0.270 0.340 0.001 0.089 1.107
beta1_yellow[12] 0.000 0.000 0.000 0.000 0.000
beta1_yellow[13] 1.740 1.517 0.068 1.102 5.703
beta1_yellow[14] 0.737 1.046 0.012 0.439 4.319
beta1_yellow[15] 0.186 0.245 0.011 0.082 0.923
beta1_yellow[16] 0.000 0.000 0.000 0.000 0.000
beta2_yellow[1] 0.000 0.000 0.000 0.000 0.000
beta2_yellow[2] -7.622 5.816 -21.956 -6.300 -0.542
beta2_yellow[3] -7.581 5.823 -22.272 -6.216 -0.403
beta2_yellow[4] -5.598 5.918 -20.326 -3.620 -0.124
beta2_yellow[5] 0.000 0.000 0.000 0.000 0.000
beta2_yellow[6] 0.000 0.000 0.000 0.000 0.000
beta2_yellow[7] 0.000 0.000 0.000 0.000 0.000
beta2_yellow[8] 0.000 0.000 0.000 0.000 0.000
beta2_yellow[9] 0.000 0.000 0.000 0.000 0.000
beta2_yellow[10] 0.000 0.000 0.000 0.000 0.000
beta2_yellow[11] -6.842 5.237 -20.016 -5.681 -0.411
beta2_yellow[12] 0.000 0.000 0.000 0.000 0.000
beta2_yellow[13] -5.189 5.301 -18.787 -3.031 -0.277
beta2_yellow[14] -6.598 4.883 -18.747 -5.884 -0.436
beta2_yellow[15] -7.500 5.586 -20.823 -6.429 -0.252
beta2_yellow[16] 0.000 0.000 0.000 0.000 0.000
beta3_yellow[1] 25.000 0.000 25.000 25.000 25.000
beta3_yellow[2] 28.853 1.090 27.305 28.661 31.696
beta3_yellow[3] 32.396 3.206 23.371 32.599 38.267
beta3_yellow[4] 30.104 3.289 23.940 29.394 35.691
beta3_yellow[5] 25.000 0.000 25.000 25.000 25.000
beta3_yellow[6] 25.000 0.000 25.000 25.000 25.000
beta3_yellow[7] 25.000 0.000 25.000 25.000 25.000
beta3_yellow[8] 25.000 0.000 25.000 25.000 25.000
beta3_yellow[9] 25.000 0.000 25.000 25.000 25.000
beta3_yellow[10] 25.000 0.000 25.000 25.000 25.000
beta3_yellow[11] 35.198 3.741 30.306 34.173 43.347
beta3_yellow[12] 25.000 0.000 25.000 25.000 25.000
beta3_yellow[13] 32.476 2.545 29.112 32.276 37.578
beta3_yellow[14] 33.832 3.432 29.187 33.497 42.516
beta3_yellow[15] 36.409 3.941 30.243 36.111 43.576
beta3_yellow[16] 35.000 0.000 35.000 35.000 35.000
beta0_black[1] 0.011 0.075 -0.137 0.011 0.159
beta0_black[2] 1.855 0.032 1.791 1.855 1.917
beta0_black[3] 1.239 0.037 1.165 1.240 1.311
beta0_black[4] 2.423 0.048 2.330 2.422 2.520
beta0_black[5] 5.137 2.370 1.875 4.675 11.386
beta0_black[6] 5.259 2.237 2.570 4.700 11.359
beta0_black[7] 4.102 2.094 1.651 3.529 9.949
beta0_black[8] 0.899 0.226 0.484 0.888 1.380
beta0_black[9] 2.808 0.262 2.314 2.803 3.318
beta0_black[10] 1.286 0.050 1.189 1.285 1.386
beta0_black[11] 2.978 0.052 2.876 2.978 3.079
beta0_black[12] 4.913 0.192 4.539 4.912 5.297
beta0_black[13] 0.458 0.337 0.017 0.390 1.228
beta0_black[14] 2.850 0.108 2.642 2.850 3.064
beta0_black[15] 0.721 0.052 0.622 0.721 0.822
beta0_black[16] 3.772 0.073 3.628 3.773 3.915
beta2_black[1] 0.000 0.000 0.000 0.000 0.000
beta2_black[2] 0.000 0.000 0.000 0.000 0.000
beta2_black[3] 0.000 0.000 0.000 0.000 0.000
beta2_black[4] 0.000 0.000 0.000 0.000 0.000
beta2_black[5] 0.000 0.000 0.000 0.000 0.000
beta2_black[6] 0.000 0.000 0.000 0.000 0.000
beta2_black[7] 0.000 0.000 0.000 0.000 0.000
beta2_black[8] 0.000 0.000 0.000 0.000 0.000
beta2_black[9] 0.000 0.000 0.000 0.000 0.000
beta2_black[10] 0.000 0.000 0.000 0.000 0.000
beta2_black[11] 0.000 0.000 0.000 0.000 0.000
beta2_black[12] 0.000 0.000 0.000 0.000 0.000
beta2_black[13] -2.621 1.809 -6.963 -2.202 -0.266
beta2_black[14] 0.000 0.000 0.000 0.000 0.000
beta2_black[15] 0.000 0.000 0.000 0.000 0.000
beta2_black[16] 0.000 0.000 0.000 0.000 0.000
beta3_black[1] 25.000 0.000 25.000 25.000 25.000
beta3_black[2] 25.000 0.000 25.000 25.000 25.000
beta3_black[3] 25.000 0.000 25.000 25.000 25.000
beta3_black[4] 25.000 0.000 25.000 25.000 25.000
beta3_black[5] 25.000 0.000 25.000 25.000 25.000
beta3_black[6] 25.000 0.000 25.000 25.000 25.000
beta3_black[7] 25.000 0.000 25.000 25.000 25.000
beta3_black[8] 25.000 0.000 25.000 25.000 25.000
beta3_black[9] 25.000 0.000 25.000 25.000 25.000
beta3_black[10] 25.000 0.000 25.000 25.000 25.000
beta3_black[11] 25.000 0.000 25.000 25.000 25.000
beta3_black[12] 25.000 0.000 25.000 25.000 25.000
beta3_black[13] 37.337 3.315 27.221 38.112 40.912
beta3_black[14] 25.000 0.000 25.000 25.000 25.000
beta3_black[15] 25.000 0.000 25.000 25.000 25.000
beta3_black[16] 25.000 0.000 25.000 25.000 25.000
beta4_black[1] -0.328 0.092 -0.508 -0.327 -0.146
beta4_black[2] 0.257 0.053 0.155 0.257 0.363
beta4_black[3] -0.720 0.064 -0.842 -0.721 -0.598
beta4_black[4] 0.636 0.129 0.388 0.634 0.898
beta4_black[5] 0.759 1.585 -1.141 0.455 4.313
beta4_black[6] 0.850 1.474 -0.884 0.533 4.527
beta4_black[7] 0.768 1.460 -0.951 0.490 3.983
beta4_black[8] 0.078 0.297 -0.564 0.091 0.655
beta4_black[9] 1.031 0.821 -0.054 0.858 3.183
beta4_black[10] 0.076 0.066 -0.057 0.077 0.203
beta4_black[11] -0.609 0.066 -0.738 -0.609 -0.481
beta4_black[12] 1.243 0.300 0.679 1.237 1.856
beta4_black[13] -1.005 0.077 -1.164 -1.004 -0.858
beta4_black[14] 0.222 0.184 -0.138 0.221 0.592
beta4_black[15] 0.003 0.047 -0.089 0.003 0.093
beta4_black[16] -0.416 0.113 -0.634 -0.417 -0.192
beta0_dsr[11] -0.329 0.009 -0.347 -0.329 -0.311
beta0_dsr[12] 1.408 0.043 1.326 1.406 1.495
beta0_dsr[13] 0.155 0.021 0.113 0.155 0.198
beta0_dsr[14] -0.127 0.017 -0.159 -0.127 -0.094
beta0_dsr[15] 0.205 0.014 0.179 0.205 0.232
beta0_dsr[16] 0.157 0.016 0.127 0.157 0.188
beta1_dsr[11] 0.000 0.000 0.000 0.000 0.000
beta1_dsr[12] 0.000 0.000 0.000 0.000 0.000
beta1_dsr[13] 0.000 0.000 0.000 0.000 0.000
beta1_dsr[14] 0.000 0.000 0.000 0.000 0.000
beta1_dsr[15] 0.000 0.000 0.000 0.000 0.000
beta1_dsr[16] 0.000 0.000 0.000 0.000 0.000
beta4_dsr[11] 0.543 0.035 0.473 0.543 0.613
beta4_dsr[12] 0.808 0.226 0.385 0.803 1.268
beta4_dsr[13] -0.078 0.041 -0.159 -0.077 0.000
beta4_dsr[14] -0.208 0.071 -0.348 -0.208 -0.067
beta4_dsr[15] 0.217 0.022 0.175 0.218 0.260
beta4_dsr[16] -0.047 0.045 -0.136 -0.047 0.041
beta0_slope[11] -2.209 0.021 -2.251 -2.208 -2.168
beta0_slope[12] -4.985 0.188 -5.379 -4.980 -4.626
beta0_slope[13] -1.501 0.034 -1.567 -1.501 -1.437
beta0_slope[14] -2.712 0.053 -2.822 -2.710 -2.613
beta0_slope[15] -1.592 0.022 -1.635 -1.592 -1.551
beta0_slope[16] -2.653 0.042 -2.734 -2.655 -2.569
beta1_slope[11] 0.000 0.000 0.000 0.000 0.000
beta1_slope[12] 0.000 0.000 0.000 0.000 0.000
beta1_slope[13] 0.000 0.000 0.000 0.000 0.000
beta1_slope[14] 0.000 0.000 0.000 0.000 0.000
beta1_slope[15] 0.000 0.000 0.000 0.000 0.000
beta1_slope[16] 0.000 0.000 0.000 0.000 0.000
beta4_slope[11] -0.335 0.075 -0.490 -0.335 -0.191
beta4_slope[12] -4.667 0.793 -6.327 -4.606 -3.266
beta4_slope[13] 0.328 0.054 0.225 0.328 0.431
beta4_slope[14] -0.158 0.152 -0.456 -0.156 0.137
beta4_slope[15] -0.360 0.033 -0.424 -0.360 -0.294
beta4_slope[16] -0.139 0.102 -0.344 -0.136 0.053
beta5_slope[11] 3.069 0.107 2.868 3.063 3.291
beta5_slope[12] 2.008 0.369 1.325 1.995 2.762
beta5_slope[13] 0.881 0.079 0.729 0.879 1.038
beta5_slope[14] 3.761 0.365 3.099 3.740 4.543
beta5_slope[15] 1.310 0.073 1.168 1.308 1.456
beta5_slope[16] 2.731 0.230 2.304 2.716 3.206
tau_comp[1] 3.780 0.295 3.221 3.778 4.381
tau_comp[2] 2.859 0.257 2.382 2.852 3.403
tau_comp[3] 0.869 0.074 0.733 0.866 1.020
tau_comp[4] 1.779 0.274 1.304 1.758 2.377
tau_comp[5] 1.262 0.183 0.941 1.249 1.654
sigma_H[1] 0.205 0.048 0.121 0.203 0.308
sigma_H[2] 0.175 0.029 0.127 0.172 0.237
sigma_H[3] 0.188 0.040 0.115 0.186 0.272
sigma_H[4] 0.462 0.077 0.330 0.454 0.636
sigma_H[5] 0.956 0.206 0.582 0.947 1.399
sigma_H[6] 0.510 0.188 0.139 0.501 0.907
sigma_H[7] 0.312 0.062 0.214 0.305 0.454
sigma_H[8] 0.425 0.087 0.291 0.415 0.621
sigma_H[9] 0.516 0.117 0.328 0.502 0.778
sigma_H[10] 0.209 0.039 0.144 0.206 0.294
sigma_H[11] 0.268 0.044 0.193 0.263 0.364
sigma_H[12] 0.381 0.144 0.197 0.337 0.725
sigma_H[13] 0.218 0.037 0.156 0.215 0.299
sigma_H[14] 0.499 0.088 0.346 0.493 0.689
sigma_H[15] 0.231 0.038 0.168 0.226 0.318
sigma_H[16] 0.216 0.040 0.148 0.212 0.306
lambda_H[1] 2.186 3.740 0.080 0.998 11.214
lambda_H[2] 5.024 4.537 0.504 3.743 17.014
lambda_H[3] 5.305 8.178 0.246 2.695 26.661
lambda_H[4] 0.006 0.004 0.001 0.005 0.017
lambda_H[5] 3.367 10.883 0.027 0.650 24.778
lambda_H[6] 3.837 9.178 0.008 0.117 32.088
lambda_H[7] 0.012 0.008 0.002 0.010 0.035
lambda_H[8] 7.521 10.109 0.077 3.933 36.315
lambda_H[9] 0.015 0.010 0.003 0.013 0.041
lambda_H[10] 0.229 0.330 0.025 0.153 0.843
lambda_H[11] 0.237 0.362 0.011 0.106 1.167
lambda_H[12] 3.830 5.103 0.184 2.207 17.023
lambda_H[13] 3.076 2.762 0.254 2.328 10.405
lambda_H[14] 3.320 4.185 0.214 2.048 13.862
lambda_H[15] 0.019 0.047 0.003 0.014 0.057
lambda_H[16] 0.586 1.034 0.035 0.282 3.052
mu_lambda_H[1] 4.052 1.888 1.038 3.844 8.116
mu_lambda_H[2] 3.499 1.928 0.440 3.282 7.811
mu_lambda_H[3] 3.394 1.810 0.739 3.140 7.545
sigma_lambda_H[1] 8.032 4.321 1.743 7.277 17.824
sigma_lambda_H[2] 7.789 4.734 0.760 6.959 18.246
sigma_lambda_H[3] 6.174 3.903 0.982 5.399 15.906
beta_H[1,1] 6.501 1.272 3.378 6.708 8.454
beta_H[2,1] 9.778 0.602 8.497 9.814 10.901
beta_H[3,1] 7.819 0.824 5.779 7.937 9.161
beta_H[4,1] 10.442 7.486 -4.289 10.427 24.719
beta_H[5,1] -0.319 2.561 -5.577 -0.165 4.194
beta_H[6,1] 2.539 4.490 -7.827 3.836 8.280
beta_H[7,1] 1.142 5.676 -10.980 1.550 11.348
beta_H[8,1] 1.474 4.662 -2.705 1.233 3.865
beta_H[9,1] 13.665 5.655 2.708 13.545 24.920
beta_H[10,1] 6.865 1.866 2.977 6.928 10.557
beta_H[11,1] 4.911 3.537 -2.671 5.620 9.994
beta_H[12,1] 2.596 1.082 0.659 2.524 4.934
beta_H[13,1] 9.105 0.891 7.220 9.167 10.589
beta_H[14,1] 2.141 1.028 0.134 2.146 4.147
beta_H[15,1] -7.108 3.356 -13.199 -7.305 -0.095
beta_H[16,1] 3.839 2.900 -0.723 3.401 10.311
beta_H[1,2] 7.778 0.268 7.232 7.785 8.283
beta_H[2,2] 9.880 0.165 9.556 9.879 10.197
beta_H[3,2] 8.799 0.207 8.385 8.803 9.216
beta_H[4,2] 3.030 1.419 0.289 3.018 5.817
beta_H[5,2] 1.736 0.991 -0.194 1.736 3.655
beta_H[6,2] 5.854 1.136 3.350 6.009 7.682
beta_H[7,2] 2.418 1.059 0.464 2.364 4.619
beta_H[8,2] 2.956 1.277 1.063 3.162 4.292
beta_H[9,2] 3.171 1.078 1.111 3.137 5.371
beta_H[10,2] 8.025 0.371 7.250 8.035 8.711
beta_H[11,2] 9.821 0.636 8.867 9.706 11.210
beta_H[12,2] 3.941 0.375 3.262 3.924 4.718
beta_H[13,2] 9.157 0.246 8.702 9.148 9.654
beta_H[14,2] 4.028 0.339 3.389 4.024 4.719
beta_H[15,2] 11.566 0.606 10.307 11.579 12.727
beta_H[16,2] 4.492 0.823 2.998 4.462 6.213
beta_H[1,3] 8.462 0.262 7.985 8.455 8.989
beta_H[2,3] 9.854 0.150 9.562 9.853 10.149
beta_H[3,3] 9.459 0.168 9.141 9.453 9.808
beta_H[4,3] -2.199 0.855 -3.909 -2.194 -0.580
beta_H[5,3] 3.906 0.595 2.679 3.911 5.071
beta_H[6,3] 8.543 1.201 6.710 8.375 10.916
beta_H[7,3] -2.572 0.634 -3.804 -2.563 -1.336
beta_H[8,3] 5.329 0.587 4.695 5.229 6.610
beta_H[9,3] -2.696 0.697 -4.121 -2.677 -1.401
beta_H[10,3] 8.560 0.278 8.016 8.564 9.111
beta_H[11,3] 8.502 0.293 7.895 8.521 9.005
beta_H[12,3] 5.272 0.321 4.528 5.313 5.776
beta_H[13,3] 8.849 0.176 8.491 8.854 9.180
beta_H[14,3] 5.785 0.267 5.213 5.807 6.248
beta_H[15,3] 10.277 0.301 9.696 10.274 10.887
beta_H[16,3] 5.986 0.614 4.741 6.001 7.120
beta_H[1,4] 8.096 0.187 7.710 8.109 8.432
beta_H[2,4] 9.844 0.130 9.585 9.850 10.080
beta_H[3,4] 9.940 0.156 9.594 9.956 10.205
beta_H[4,4] 11.815 0.468 10.884 11.809 12.715
beta_H[5,4] 5.596 0.751 4.360 5.519 7.285
beta_H[6,4] 7.027 0.943 5.191 7.114 8.501
beta_H[7,4] 8.547 0.354 7.874 8.545 9.258
beta_H[8,4] 6.707 0.282 6.126 6.731 7.136
beta_H[9,4] 7.423 0.442 6.554 7.420 8.286
beta_H[10,4] 7.519 0.224 7.089 7.512 7.992
beta_H[11,4] 9.504 0.195 9.132 9.505 9.878
beta_H[12,4] 7.254 0.198 6.880 7.251 7.667
beta_H[13,4] 9.063 0.141 8.781 9.068 9.333
beta_H[14,4] 7.838 0.215 7.436 7.834 8.276
beta_H[15,4] 9.482 0.225 9.046 9.482 9.922
beta_H[16,4] 9.403 0.246 8.940 9.397 9.902
beta_H[1,5] 8.899 0.151 8.579 8.904 9.195
beta_H[2,5] 10.674 0.101 10.484 10.671 10.882
beta_H[3,5] 10.765 0.166 10.470 10.751 11.101
beta_H[4,5] 8.331 0.490 7.369 8.320 9.334
beta_H[5,5] 5.174 0.613 3.754 5.252 6.154
beta_H[6,5] 9.034 0.621 8.048 8.947 10.335
beta_H[7,5] 6.576 0.348 5.896 6.574 7.245
beta_H[8,5] 8.248 0.229 7.886 8.229 8.724
beta_H[9,5] 8.206 0.454 7.312 8.201 9.091
beta_H[10,5] 10.046 0.231 9.591 10.047 10.479
beta_H[11,5] 11.590 0.228 11.141 11.594 12.034
beta_H[12,5] 8.614 0.188 8.240 8.611 8.987
beta_H[13,5] 10.140 0.134 9.883 10.135 10.416
beta_H[14,5] 9.252 0.218 8.844 9.244 9.701
beta_H[15,5] 11.310 0.234 10.843 11.312 11.787
beta_H[16,5] 9.839 0.187 9.450 9.842 10.198
beta_H[1,6] 10.118 0.236 9.737 10.086 10.688
beta_H[2,6] 11.475 0.118 11.236 11.475 11.700
beta_H[3,6] 10.612 0.159 10.277 10.620 10.901
beta_H[4,6] 12.917 0.884 11.135 12.955 14.667
beta_H[5,6] 5.684 0.624 4.579 5.659 7.007
beta_H[6,6] 8.707 0.585 7.435 8.738 9.783
beta_H[7,6] 10.552 0.603 9.403 10.532 11.721
beta_H[8,6] 9.510 0.302 8.984 9.533 9.972
beta_H[9,6] 8.261 0.773 6.764 8.258 9.830
beta_H[10,6] 9.136 0.344 8.376 9.162 9.734
beta_H[11,6] 10.688 0.367 9.894 10.712 11.344
beta_H[12,6] 9.385 0.238 8.967 9.372 9.902
beta_H[13,6] 11.086 0.156 10.781 11.082 11.402
beta_H[14,6] 9.923 0.286 9.331 9.931 10.461
beta_H[15,6] 10.553 0.416 9.726 10.555 11.370
beta_H[16,6] 10.600 0.258 10.052 10.603 11.091
beta_H[1,7] 10.162 1.266 6.833 10.437 11.939
beta_H[2,7] 12.106 0.524 10.988 12.139 13.062
beta_H[3,7] 10.283 0.688 8.783 10.328 11.480
beta_H[4,7] 1.398 4.534 -7.517 1.408 10.208
beta_H[5,7] 5.890 2.036 1.955 5.922 9.889
beta_H[6,7] 6.743 3.355 -1.307 7.741 11.084
beta_H[7,7] 7.674 3.047 1.812 7.697 13.510
beta_H[8,7] 10.927 1.154 9.170 10.873 12.805
beta_H[9,7] 5.018 4.006 -3.178 5.022 12.588
beta_H[10,7] 10.079 1.612 7.314 9.915 13.728
beta_H[11,7] 11.370 1.776 8.262 11.207 15.352
beta_H[12,7] 9.795 0.951 7.707 9.884 11.428
beta_H[13,7] 11.703 0.712 10.044 11.770 12.909
beta_H[14,7] 10.632 0.941 8.656 10.671 12.385
beta_H[15,7] 13.961 2.192 9.607 13.974 18.213
beta_H[16,7] 12.781 1.382 10.471 12.655 15.910
beta0_H[1] 8.879 17.454 -27.803 8.891 47.069
beta0_H[2] 10.382 7.727 -5.007 10.394 26.929
beta0_H[3] 9.256 10.151 -12.403 9.436 29.919
beta0_H[4] 0.734 186.498 -361.080 -0.853 371.596
beta0_H[5] 4.105 29.020 -50.014 4.388 56.644
beta0_H[6] 7.992 57.690 -123.755 7.799 133.766
beta0_H[7] 6.341 135.580 -270.503 5.507 279.325
beta0_H[8] 7.678 38.197 -18.926 6.531 35.883
beta0_H[9] 6.216 122.080 -237.789 6.201 251.139
beta0_H[10] 7.885 38.485 -71.999 8.063 84.549
beta0_H[11] 10.622 53.979 -98.936 10.597 123.651
beta0_H[12] 6.666 12.963 -18.108 6.545 31.515
beta0_H[13] 9.704 11.287 -12.093 9.779 30.364
beta0_H[14] 7.331 11.413 -15.018 7.142 30.301
beta0_H[15] 10.111 112.798 -224.869 11.319 227.900
beta0_H[16] 9.284 29.495 -53.377 8.952 72.667